Phononic crystals with planar defects 
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We study the effect of planar defects in phononic crys- 
tals of spherical scatterers. It is shown that a plane of im- 
purity spheres introduces modes of vibration of the elastic 
field localized on this plane at frequencies within a fre- 
quency gap of a pure phononic crystal; these show up as 
sharp resonances in the transmittance of elastic waves in- 
cident on a slab of the crystal. A periodic arrangement of 
impurity planes along a given direction creates narrow im- 
purity bands with a width which depends on the position 
of these bands within the frequency gap of the pure crystal 
and on the separation between the impurity planes. We 
show how a slight deviation from periodicity (one impu- 
rity plane is different from the rest) reduces dramatically 
the transmittance of elastic waves incident on a slab of the 
crystal. 
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I. INTRODUCTION 

Phononic crystals are composite materials whose elas- 
tic properties, as described by the mass density p and 
the Lame coefficients A and /i, vary periodically in space. 
In recent years there has been a growing interest in the 
study of phononic crystals, especially in relation to the 
possibility of phononic band gap materials (see, e.g., 
Ref. |l| and references therein). By definition these ex- 
hibit regions of frequency (phononic gaps) over which no 
vibrations are possible within the infinite crystal. This 
implies that an elastic wave, whatever its direction of 
propagation and its polarization (longitudinal or trans- 
verse), incident on a slab of the material of some thickness 
will be totally reflected by it, if its frequency lies within 
a phononic gap. Obvious technological applications of 
the above are non-absorbing mirrors and vibration-free 
cavities which might be very useful in high-precision me- 
chanical systems operating in a given frequency range. 

On the theory side one would like to be able to calcu- 
late the frequency band structure of a phononic crystal 
and also the transmission and reflection coefficients of 
elastic waves incident on a slab of rthe material of finite 
thickness. In a previous publication^ (to be referred to as 
paper I) we presented a formalism which can do that for 
phononic crystals consisting of non-overlapping spheres 
embedded in an elastic medium of different mass density 
and Lame coefficients. In the present paper we apply the 



formalism of paper I to a specific system consisting of 
non-overlapping solid spheres embedded in a solid host 
medium. We assume that the spheres are centered on the 
sites of an fee lattice and view the crystal as a sequence of 
(001) planes along the z axis. We calculate the frequency 
band structure of the infinite crystal, and the transmis- 
sion coefficient for an elastic wave incident on a slab of 
the material parallel to the (GDI) surface consisting of a 
finite number of planes of spheres. For a thick slab the 
transmittance vanishes over the exact range of the band 
gap (of the infinite crystal) as one would expect; but of 
greater interest is the fact that the transmittance behaves 
in almost the same manner for a thin slab consisting of 
just two planes of spheres. 

Linear and point defects in phononic crystals have been 
investigated by Torres et al.n In the present work we 
study the effect of planar defects, i.e. of impurity planes, 
on the transmittance of a slab. An impurity plane has 
the same two-dimensional (2D) periodicity as the rest of 
the planes of the slab, but the spheres centered on this 
plane are different. In our example they have a different 
radius. This leads to gap states: vibrational modes of the 
elastic field localized on the impurity plane, which show 
up as transmission resonances at frequencies within the 
gap of the pure crystal. 

Finally we consider a crystal whose unit cell along the z 
direction consists of a number of planes one of which (im- 
purity plane) is different from the other planes of the unit 
cell. This leads to impurity bands within the frequency 
gap of the pure crystal, the width of which depends on the 
position of the impurity band within the gap and on the 
separation between the impurity planes. We show that 
a slight deviation from periodicity reduces dramatically 
the transmission coefficient of an elastic wave incident on 
a slab of the crystal. 

In all calculations of the transmittance of a slab we 
assume that the material that exists between the spheres 
of the slab extends to infinity on either side of the slab. 



II. PROPERTIES OF THE PURE CRYSTAL 

The pure crystal consists of non-overlapping lead 
spheres centered on the sites of an fee lattice. The elastic 
properties of the spheres are characterized by the mass 
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density, p, and the longitudinal and transverse velocities 
of propagation, denoted by c; and ct respectively. The 
latter are given by: = (A + 2p)/p and Cj — ^1 p. For 
lead we have: p ~ 11.357 gr cm^^, c; — 2158 ms~^, 
Ct — 860 ms^^. The space between the spheres is 
filled with epoxy ( p = 1.180 gr cm^'^, q — 2540 m s^^, 
Ct = 1160 m s-i). 

We calculate the complex frequency band structure of 
the infinite crystal (extending from z — oo to z — > oo) 
associated with the (001) surface in the manner described 
in paper I. For any given reduced wavevector k|| , within 
the surface Brillouin zone (SBZ) of this surface, we calcu- 
late the frequency lines kz(uj] k||) as functions of the an- 
gular frequency to. Our results, obtained with an angular 
momentum cut-off imax = 4 and 13 2D reciprocal vec- 
tors, g, are converged within an accuracy 10~^. The or- 
dinary frequency band structure which interests us here, 
corresponds to real sections (regions of frequency over 
which fcz([j;k||) is real) of these lines. In Fig. |l](a) we 
show the frequency bands of the elastic field for k|| = 
(along the direction normal to the (001) plane). We 
present our results in dimensionless units; a denotes the 
lattice constant of the fee lattice and Ct the propaga- 
tion velocity of transverse elastic waves in epoxy. The 
radius S of the spheres equals 0.25a in the present case 
which corresponds to a fractional volume occupied by the 
spheres, denoted by /, equal to 0.262. 
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FIG. 1. The phononic band structure at the center 
(k|| = 0) of the SBZ of a (001) surface of an fee crystal 
of lead spheres in epoxy (a); and the corresponding trans- 
mittance curve of a slab of 16 planes of spheres parallel to 
the same surface (b). The fractional volume occupied by the 
spheres is 0.262. In (a) the black lines represent longitudinal 
modes, the grey lines transverse modes, and the dotted lines 
are deaf bands. Correspondingly, in (b) the black line shows 
the transmittance of longitudinal incident elastic waves, and 
the shaded curve that of transverse incident waves. 

The bands shown in Fig. |l|(a) by black (grey) lines 
correspond to longitudinal (transverse) modes of vibra- 
tion of the elastic field inside a sufficiently thick slab of 
the material. These modes are excited only by longitu- 
dinal (transverse) waves incident normally on the sur- 



face of the slab. This is demonstrated quite clearly in 
Fig. |l|(b), which shows the transmission coefficient of a 
plane wave incident normally on a slab of the material 
consisting of 16 planes of spheres parallel to the (001) 
surface. The shaded curve shows the transmission coef- 
ficient for a transverse incident wave and the black line 
that of a longitudinal incident wave. 
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FIG. 2. Transmittance of longitudinal (black lines) and 
transverse (shaded curves) elastic waves incident normally on 
a slab of the phononic crystal described in the caption of 
Fig. |l| consisting of two planes of spheres (a), and 16 planes 
of spheres (b). 

We hasten to add (see also paper I) that in a com- 
posite material the modes of vibration of the elastic field 
need not be purely longitudinal or purely transverse even 
along symmetry directions and, therefore, the fact that 
they behave so in the system under consideration is worth 
noting. Off the normal direction (k|| ^ 0) the frequency 
bands are clearly hybridized and they are excited by ei- 
ther longitudinal or transverse incident plane waves. The 
bands shown by dotted lines in Fig. |l|(a) are deaf bands. 
The g = component of the corresponding modes of 
vibration vanishes and that being the only component 
which couples with the external elastic field in the fre- 
quency range under consideration, they cannot be excited 
by an incident plane wave, leading to total reflection of 
the latter. Finally, we point out the oscillations in the 
transmission coefficient, over the allowed regions of fre- 
quency, which are due to interference effects resulting 
from multiple reflection between the surfaces of the slab 
(Fabry-Perrot-type oscillations) . 

In Fig. 1^ we show the transmission coefficient of an 
elastic wave incident normally on a thin (001) slab con- 
sisting of just two planes of spheres, and compare it 
with the same quantity for a thick slab (of 16 planes of 
spheres) , which as we have seen is determined by the fre- 
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quency band structure of the infinite crystal. We see that 
over the frequency region that constitutes the frequency 
gap of the infinite crystal, the transmittance of a slab con- 
sisting of just two planes of spheres is not zero, but it is 
nevertheless about two orders of magnitude smaller than 
its value at frequencies below or above the gap edges. 
The Fabry-Perrot-type oscillations are of course differ- 
ent in the two cases since they depend directly on the 
thickness of the slab. 
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FIG. 3. Projection of the frequency band structure on the 
SBZ of the (001) surface of the fee phononic crystal described 
in the caption of Fig. 0. The blank areas show the frequency 
gaps in the considered frequency region. The inset shows the 
SBZ of the (001) fee surface. 

An absolute frequency gap extends between uia/ct — 
5.29 and w„a/ct = 7.38, as demonstrated in Fig. |3[ This 
shows the projection of the frequency band structure on 
the SBZ, along the symmetry lines of the latter shown in 
the inset. The shaded areas correspond to frequencies for 
which there exists at least one propagating Bloch wave 
(an eigenmode of the elastic field) in the infinite crys- 
tal; in the frequency regions represented by the blank 
areas none exists. We verified, by calculating the fre- 
quency lines at selected points of the SBZ, that the gap 
does exist over the entire SBZ, and that it is therefore an 
absolute phononic gap extending from w/a/ct = 5.29 to 
LOua/ct = 7.38. The above values of the lower and upper 
edges of the gap are different from the values ujia/ct ~ 6.2 
and] ujua/ct ~ 8.0, respectively, found by Kafesaki et 
al.p who studied the same system using the plane- wave 
method. We believe our results to be the more accu- 
rate. It is now recognized that the plane-wave method is 
not particularly suited to the study of phononic crystals 
consisting of non-overlapping spheres in a host medium, 
because of the very large nuiaber of plane waves required 
to obtain convergent results. □ 

III. PLANAR DEFECTS 

Next we consider a slab of the phononic crystal consist- 
ing of five planes of spheres, one of which may be different 



from the other four. It has the same 2D periodicity par- 
allel to the surface of the slab, but the spheres of this, 
so-called, impurity plane may be different: they have a 
smaller or larger radius than the spheres of the other 
planes. Our results are summarized in the six diagrams 
of Fig. ^. The first three, (a), (b) and (c), show the trans- 
mittance of a transverse elastic wave incident normally 
on a slab of five planes of spheres. For Fig. ^(a) all five 
planes are the same, as in the infinite crystal described 
in the preceding paragraphs, and the transmittance van- 
ishes, as expected, over the frequency region correspond- 
ing to the frequency gap of the infinite crystal for k|| — 0, 
shown in Fig. |l|(a). For Figs. ||(b) and |(c) the middle 
plane is different from the other four. The spheres of 
this plane have a smaller radius: Si — 0.8S for Fig. ^(b) 
and Si = 0.65* for Fig. ^(c), where S is the radius of the 
spheres in the other planes. The transmission resonance 
that now appears at a frequency within the gap signifies 
the existence of a state of the elastic field centered on 
the impurity plane: a mode of vibration of the elastic 
field that extends to infinity parallel to the surface of the 
slab (in the manner of a Bloch wave), but decays rapidly 
normal to the impurity plane on either side of it. It ap- 
pears that the normal mode of vibration (with ky = 0) at 
the top of the valence band (we use the term by analogy 
to semiconductor physics to denote the frequency band 
below the gap), splits off this band, becoming a local- 
ized (on the impurity plane) vibration, with a frequency 
higher in the gap the smaller the radius of the impurity 
spheres. Figs, ^(d) to ^(f) demonstrate the same phe- 
nomenon for longitudinal waves incident normally on the 
slab. For Fig. ^(d) all five planes are the same as in the 
infinite crystal, and the transmission coefficient accords 
with the corresponding frequency band structure shown 
in Fig. [|(a). For Figs. |(e) and |(f) the spheres of the 
middle plane have a radius smaller than that of the other 
four planes: Si = 0.6S for Fig. |(e) and S", = 0.55" for 
Fig. H(f). We note that while the trend is the same, the 
resonance frequency stays closer to the top of the valence 
band for longitudinal waves. 

We note that the transmission coefficient equals unity 
at the resonance frequency. This is the case when the 
impurity plane lies in the middle of the slab; when the 
impurity plane is any other than the middle plane, the 
transmission coefficient at the resonance frequency is less 
than unity, as demonstrated in Fig. |5| This figure shows 
the transmission coefficient for a transverse plane wave 
incident normally on a slab of five planes of spheres one 
of which is an impurity plane. The spheres centered on 
this plane have a radius Si = 0.75", where S is the ra- 
dius of the spheres of the other planes. The top dia- 
gram is obtained when the impurity plane is the middle 
plane of the slab, the middle diagram is obtained when 
the impurity plane is the second from the surface, and 
the third diagram is obtained when the impurity plane 
lies at the surface of the slab. Only in the first case is 
the transmission coefficient at resonance equal to unity. 
When the impurity plane is removed from the center of 
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the slab by one plane, the value of the transmission coef- 
ficient at resonance diminishes by at least two orders of 
magnitude (middle diagram), and the resonance disap- 
pears altogether when the impurity plane is removed to 
the surface of the slab. Similar results are obtained for a 
longitudinal incident wave. This is indeed a general char- 
acteristic of resonant tunneling. We note, for example, 
the similarity of the above results with those obtained in 
the transmission of an electron through a double barrier, 
when the well between the two barrier tops allows for 
a resonance state of the tunneling electron. When the 
double barrier is symmetric the transmission coefRcient 
at resonance equals unity, when the double barrier is not 
symmetric the transmission coefficient is less than unity.El 
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waves ( the displacement vector is parallel to the surface 
of the slab); the third band is excited by both longitudi- 
nal and p-polarized waves (the displacement vector lies 
in the plane of incidence). We note that all three bands 
of gap states are restricted within a rather narrow range 
of frequencies about the midgap frequency. 
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FIG. 5. Transmittance of transverse elastic waves incident 
normally on a slab of an fee crystal of lead spheres in epoxy. 
The slab consists of five planes of spheres parallel to the (001) 
surface. The spheres have a radius S = 0.25a, except those of 
the middle plane (top diagram), or of the second plane from 
the surface (middle diagram), or of the surface of the slab 
(bottom diagram), which have Si — 0.7S . 



FIG. 4. Transmittance of elastic waves incident normally 
on a slab of an fee crystal of lead spheres in epoxy. The slab 
consists of five planes of spheres parallel to the (001) surface. 
The spheres have a radius S = 0.25a, except those of the 
middle plane which have a different radius Si [(a): Si — S, 
(b): Si — 0.8S, (c): 5"; — 0.6S, for a transverse incident wave; 
(d): Si = S, (e): S^ = 0.65, (f): S^ = 0.55*, for a longitudinal 
incident wave]. 

In Fig. ^ we show the variation with k|| of the resonance 
frequency, associated with an impurity plane at the mid- 
dle of a slab of five planes of spheres. The spheres of the 
middle plane have a radius St = 0.7S, where S = 0.25a 
is the radius of the spheres in the other planes. The gap 
states/resonances along XM and Mf are hybridized, ex- 
cited by both longitudinal and transverse incident plane 
waves. The states/resonances along f X corresponding to 
the two top bands are excited by s-polarized transverse 
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FIG. 6. Variation with ky of the resonance frequency as- 
sociated with an impurity plane at the center of a slab of 
five (001) fee planes of lead spheres in epoxy {S — 0.25a, 
S^ = 0.75). 



4 



Next we consider an infinite crystal (extending from 
z — !■ — oo to z — > oo) the unit cell of which, along the 
z axis, consists of three fee planes of spheres; we have: 
• • • SSS^SSS.SSS, The spheres of the third plane 
of the unit slice have a radius Si which is smaller than 
the radius S of the spheres of the other two planes. We 
put Si = 0.7S. The resonances shown in Fig. ^ now 
develop into bands shown by shaded regions in Fig. |^. 
These bands are nevertheless extremely narrow for most 
k|| , which implies that the vibrations are strongly local- 
ized on the impurity planes, with very little interaction 
between impurity planes. This is, of course, consistent 
with the sharpness of the resonance seen in a slab of five 
planes (see top diagram of Fig. ||). 
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FIG. 7. Projection of the frequency band structure on the 
SBZ of the (001) fee surface of a phononic crystal of lead 
spheres in epoxy. The spheres are centered on the sites of an 
fee lattice and the unit cell of the structure along the z axis 
consists of three planes of spheres: SSSi, with S = 0.25a and 
Si = 0.75. The dotted lines reproduce the results shown in 
Fig. ^ for comparison. 

We can vary the width of the impurity bands by push- 
ing them away from the middle of the gap. Putting 
Si = 0.8 S displaces the impurity bands of Fig. |^ to- 
wards the bottom of the gap; at the F point this leads 
to a band- width of Auja/ct ~ 0.14 which is more than 
three times that shown in Fig. ^ which equals 0.04. We 
can understand this by noting that the interaction be- 
tween states localized on neighboring impurity planes 
depends on the separation, d, between the planes ap- 
proximately as ~ exp [— I fci (w; k|| ) I d ] , where ki (w; ky ) is 
the smallest (in magnitude) of the imaginary parts of the 
frequency lines kz{ijJ] k||) in the pure crystal over the fre- 
quency range of the impurity band. | fci(w; k||)| is largest 
in the middle of the gap (for a given k|| ) and goes to zero 
at the edges of this gap. 

Finally, in Fig. ^ (curve (A)) we show the transmission 
coefficient in the frequency region of the impurity band 
for a transverse wave incident normally on a slab of a 
phononic crystal consisting of seven unit cells along the 
z direction, when the unit shce consists of three (001) fee 
planes: SSSt, with the spheres of the third plane having 



a radius S^ = 0.7S with S = 0.25a. 
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FIG. 8. (A): Transmittance of transverse elastic waves in- 
cident normally on a slab of seven unit slices of the phononic 
crystal described in the caption of Fig. 0, over the frequency 
range of the impurity band within the absolute gap. (B): As 
in (A) except that Si of the third unit slice is a little different 
(it equals O.GS). 

We note that the frequency gap of the pure crystal 
{Si = S) for k|| = extends from uji a/ct — 5.29 to 
ujua/ct = 7.67. We recall that the gap states corre- 
sponding to the first six impurity planes (taken individ- 
ually) manifest themselves as transmission resonances at 
ujra/ct ~ 6.29 (see Fig. |^) while the seventh impurity 
plane, at the surface of the slab, does not give a trans- 
mission resonance, as shown in the bottom diagram of 
Fig. ^. The interaction between the impurity planes re- 
moves the degeneracy of the above resonance level, lead- 
ing to six discrete resonances in the transmission spec- 
trum, as shown in Fig. ||. Curve (B) of Fig. ^ is ob- 
tained for a slab which differs from the above in that the 
spheres of the third plane of the third unit slice have a 
radius Si = 0.65, instead of Si — 0.75. We note the 
dramatic reduction, by about four orders of magnitude, 
in the transmission coefficient as a result of the presence 
of just one different impurity plane. In a way this phe- 
nomenon, which may have interesting technological ap- 
plications, is to be expected; the different impurity plane 
decouples the vibrations on its right from those on its 
left. This suggests that a random distribution of impu- 
rity planes will lead to vibrational modes localized over 
smaller regions of the slab. We have here a possibility to 
study Anderson localization of classical waves in an effec- 
tively one-dimensional system (because of k|| conserva- 
tion) which allows direct comparison with corresponding 
experimental data, without the difficulty associated with 
correlation in the corresponding electronic problem. 
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